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We numerically investigate, using the time evolving block decimation algorithm, the quantum 
transport of ultra-cold bosonic atoms in a double well optical lattice through slow and periodic 
modulation of the lattice parameters (intra- and inter- well tunneling, chemical potential, etc.). The 
transport of atoms does not depend on the rate of change of the parameters (as along as the change 
is slow) and can distribute atoms in optical lattices at the quantized level without involving external 
forces. The transport of atoms depends on the atom filling in each double well and the interaction 
between atoms. In the strongly interacting region, the bosonic atoms share the same transport 
properties as non-interacting fermions with quantized transport at the half filling and no atom 
transport at the integer filling. In the weakly interacting region, the number of the transported 
atoms is proportional to the atom filling. We show the signature of the quantum transport from 
the momentum distribution of atoms that can measured in the time of flight image. A semiclassical 
transport model is developed to explain the numerically observed transport of bosonic atoms in the 
non-interacting and strongly interacting limits. The scheme may serve as an quantized battery for 
atomtronics applications. 

PACS numbers: 03.75.Lm, 03.65.Vf 



I. INTRODUCTION 

Quantum charge pumping, a coherent quantum trans- 
port process that generates steady charge currents of 
electrons through adiabatically and periodically time- 
varying potentials, is a standard method for the charge 
transport in solid state circuits I n the quantized 

charge transport, the number of particles pumped out 
during each cycle of the potential modulation is an in- 
teger and can be understood using a Berry phase [fj. 
Because of the precise control on the amount of pumped 
charges at the single electron level, the quantum charge 
pumping has found important applications in many elec- 
tronic devices @, Q . It also provides a solid foundation to 
the modern theory of electric polarization [9l-ll lj . Similar 
ideas have also been extended to another degree of free- 
dom of electrons, the spin, leading to the quantum spin 
pumping that plays a crucial role in spintronics |12l-[l4{. 
However, the exact quantized charge transport is usually 
difficult to observe in experiments because of unavoidable 
complexity due to impurities, disorder and interactions 
in the solid. 

In view of the significance of electronics and spintron- 
ics, there has been a great deal of interest recently for de- 
veloping a one-to-one analog of complex and interesting 
electronic materials, circuits, and devices using ultra-cold 
neutral atoms [15l-[l7j. This field, known as atomtronics, 
is a significant extension of the recent great efforts on em- 
ulating condensed matter physics using ultra-cold atoms 
[18|-[20j. Important concepts such as atomic batteries, 
diodes, and transistors, have been proposed recently for 
atomtronics. It would be natural and also important to 
investigate the coherent quantum transport process of 
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cold atoms in optical lattices, which is an important el- 
ement of atomtronic devices. A straightforward method 
for the transport, of course, is by applying an external 
force. However, using external forces not only involves 
the acceleration of atoms that may change the motional 
states of atoms, but also may be limited by the Bloch 
oscillation of atoms in optical lattices. Recently, a lot of 
attention has been focused on the quantum transport of 
atoms in optical lattices without involving external forces 
[2l|,[22|. For instance, it has been proposed that quantum 
transport of atoms may be achieved in optical double well 
lattices through the fast periodic modulation (in a non- 
adiabatic manner) of the inter- and intra- well tunnelings 
[2lj . It was also proposed [22j that a fast oscillating lin- 
ear potential in optical lattices can yield the quantum 
ratchet effect due to coherent destruction of tunneling 
[23|, leading to quantum transport of atoms. However, 
all these schemes require accurate control of the lattice 
parameters and their time variation. 

In this paper, we study the quantum transport of ultra- 
cold bosonic atoms in optical lattices without involving 
the accurate control of the lattice parameters and their 
time dependence. We consider atoms that are initially 
prepar ed in certain region of a double well optical lattice 
[2J-[34(. We show that quantized number of atoms can 
be transferred to another region of the lattice through 
periodic and slow modulation of the optical lattice pa- 
rameters (intra- and inter-well tunneling, chemical po- 
tential, etc.). Such quantized atom pumping may serve 
as a quantum atom battery for atomtronic devices. Note 
that there is an important difference between the atom 
transport in optical lattices and the charge transport in 
the solid. While the electrons in the solid are fermions, 
cold atoms in optical lattices can be bosons, therefore 
the Bose-Einstein statistics, instead of the Fermi-Dirac 
statistics, governs the transport dynamics. This impor- 
tant difference will be illustrated in this paper. Another 
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FIG. 1. (Color online) Schematic plot of the quantum trans- 
port of ultra-cold bosonic atoms in a double well optical lat- 
tice, (a) Illustration of the double well optical lattice pa- 
rameters, (b) Illustration of the quantum transport process 
of bosonic atoms from region A to region B in the lattice, 
(c) The periodic modulation of the lattice parameters for the 
quantum transport process. Two different types of loops in 
the parameter space are considered: square (solid line) and 
elliptical (dashed line). 



important difference is that while electrons are periodi- 
cally distributed in crystals, the density distribution of 
atoms can be prepared locally in the optical lattice and 
transferred to another region through the quantum trans- 
port process with the periodic modulation of the lattice 
parameters. We emphasize that although the periodic 
modulation of the lattice parameters in the parameter 
space are needed in the transport process, the rate of 
change of the parameters, the initial and final values of 
each parameter during its variation do not need to be 
accurately controlled. This is because the transport pro- 
cess depends on the topology of the loop in the parameter 
space, instead of the exact parameter loop or the rate of 
change of the parameters along the loop. Therefore the 
quantum transport process is very robust against errors 
in the parameter modulation. 

We find the transport of bosonic atoms depends 
strongly on the atom filling per double well as well as the 
interaction. The strongly interacting bosonic atoms be- 
have similarly as the non-interacting fermions with quan- 
tized transport at the half filling and no mass transport 
at the integer filling. In the weakly interacting region, the 
number of transported atoms is proportional to the atom 
filling. The investigation is based on the numerical sim- 
ulation of the exact quantum dynamics of cold atoms in 
double well optical lattices using the time evolving block 
decimation (TEBD) algorithm [35|, [36|. The transport 
properties observed in the numerical simulation are also 
understood by developing an analytical theory using a 
semiclassical transport model in the non-interacting and 
strongly interacting limits. 

The paper is organized as follows: Section II describes 



the physical system: the cold bosonic atoms in a double 
well optical lattice. Section III presents the numerical 
results on the transport of bosonic atoms in the optical 
lattices through the periodic modulation the lattice pa- 
rameters. In Section IV, we provide a semiclassical trans- 
port model to explain the numerically observed transport 
of bosonic cold atoms in the non-interacting and strongly 
interacting limits. Section V consists of conclusions. 



II. COLD BOSONS IN DOUBLE WELL 
OPTICAL LATTICES 



Consider ultra-cold bosonic atoms confined in a one 
dimensional double well optical lattice (Fig. la). An op- 
tical lattice is a standing wave of coherent off-resonance 
light created by the interference of two or more laser 
beams. In experiments, the ID double well optical lat- 
tice has been realized by superimposing two laser beams 
with two different wavelengths A2 = Ai/2 [24]-[27j, or the 
Fourier synthesis of asymmetric optical potentials with 
spatial periodicity \/2n (n is integer) [37]. In these ex- 
periments, the dynamics along the other two dimensions 
are frozen to the ground states by optical lattices with 
high potential depths (^> recoil energy). With the stan- 
dard tight-binding approximation, the dynamics of atoms 
in the double well optical lattices can be described by a 
Bose-Hubbard Hamiltonian [38| 
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/2 describes the inter- or intra- well 

tunneling, Vj is the external potential (e.g., the harmonic 
trap, disorder, etc.), 2 A is the chemical potential differ- 
ence between two neighboring lattice sites in a double 
well, and U is the on-site interaction strength between 

atoms. In experiments, the tunneling J — (— 1) J 5 /2 

can be adjusted by varying the intensities of the two laser 
beams, the chemical potential A can be tuned by shift- 
ing one laser beam with respect to the other, the on-site 
interaction U between atoms can be changed using the 
Feshbach resonance [39|, the disorder potential can be 
created using optical speckle potentials [ioj . 

The transport process of cold atoms is illustrated in 
Fig. lb. Initially, the bosonic atoms are confined within 
a group of lattice sites (denoted as region A) using a 
harmonic or box trapping potential [4l|, which will be 
removed after the transport starts. The transport is re- 
alized by moving the atoms in the region A to another 
group of lattice sites (denoted as region B) by modulat- 
ing the lattice parameters periodically, but without the 
actual movement of the lattice. Therefore it is a quan- 
tum transport process through the tunneling dynamics. 
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FIG. 2. (Color online) Plot of the lowest two eigenenergies 
with respect to U for two atoms in a double well. Solid line: 
A = 0. Dashed line: A = 10 J. 



To simplify the study, we assume a uniform distribution 
of the atoms in the region A at the time t = 0. During 
the transport process, the lattice parameters 5 and A are 
tuned along a periodic close loop in the parameter space, 
as illustrated in Fig. lc. 

The time-dependent dynamics of atoms in the optical 
lattice are investigated numerically using the TEBD al- 
gorithm HI, [36|. The TEBD algorithm is an effective 
method for simulating the exact one-dimensional quan- 
tum dynamics when the quantum entanglement of the 
system is low. It works well for the Hamiltonian (pQ) 
where the only non-local terms are the nearest neigh- 
boring tunneling. To check the validity of our numerical 
program, we have compared the results from the TEBD 
and that from the exact diagonalization for a small lattice 
system and they agree with a high accuracy. Henceforth, 
we set the energy unit of the system as J, and the time 
unit as J -1 . 

We assume the initial parameter 5/ J = 1 to obtain a 
well-defined localized initial state. With this parameter, 
the inter-well tunneling vanishes and the initial wave- 
function of the system is a product of the wavefunction 
in each double well. We use the number states of atoms 
at each lattice site as the basis states to represent the 
wavefunction. For instance, consider a double well lat- 
tice with a lattice length N = 20 and A/J = initially. 
Ten atoms are uniformly distributed in the region A be- 
tween sites 5 and 14. In the region A, the wavefunction 
in the /-th double well can be written as 
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where c 20 = c 02 - Wuw+u)Vuw > 

^Ju2+4j ¥~' \ nm ) ls t ne state with n (m) atoms in 
the left (right) well. Outside the region A, the wavefunc- 
tion in each double well is just <pi = |00) z . 

The quantum transport process are accomplished by 
modulating the lattice parameters (5, A) periodically in 
the parameter space, as illustrated in Fig. lc. Two differ- 
ent loops (square and elliptical) are considered and their 
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FIG. 3. (Color online) Plot of the number of atoms in the 
lattice sites 15 and 16 with respect to the atom interaction 
strength U after one cycle of the parameter modulation illus- 
trusted in the insets. Initially ten atoms uniformly occupy 
the lattice sites 5 to 14. (a) A square loop, (b) An elliptical 
loop. 



effects will be discussed later in the paper. We choose the 
rate of the change of the parameters (5, A) such that the 
tunneling process is adiabatic in each double well and 
the atoms always stay on the ground state of the dou- 
ble well. During the whole transport process, the atoms 
are assumed to be on the motional ground state at each 
lattice site. This assumption can be realized by using 
a double well lattice with an intra-well potential depth 

7 2 

Vl ~ 6Er, where Er = 2 m\ 2 1S ^ ne a ^ om rec °il energy, 
h is the Planck constant, A is the wavelength of the short 
wavelength laser, and m is the atom mass. This potential 
depth yields a tunneling rate J ~ 0.05Er (42j| . There- 
fore the required chemical potential shift A ~ 0.5 Er for 
the transport process are much smaller than the energy 
gap between the ground and first excited motional bands 
(typically several Er). 

In each isolated double well with two atoms, the Hamil- 
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FIG. 4. (Color online) Plot of the density distribution of 
atoms after 10 cycles of the parameter modulation with the 
square loop. Solid line: initial density distribution. Dashed 
line: final density distribution, (a) U — 5 J. (b) U — 20J. 
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tonian can be written as 



2A + U 






-2A + U 
-J 



(3) 



on the basis {|20) , 1 02) , |11)}. Without interaction be- 
tween atoms (U = 0), the eigenenergies are E± = 
±V4A 2 + 2 J 2 , E 3 = 0. Therefore the energy gap be- 
tween the ground state E_ and the first excited state 
Es is larger than y/2J. In the presence of strong in- 
teraction U (neglect A), the three eigenenergies are 

E± = \U (l ± y/l + 8J 2 //7 2 ) ,E 3 = U. The energy gap 

between the ground state E_ and the first excited state 
Es is ~ U. The crossover between different energy bands 
in Fig. 2 is essentially a Landau-Zener tunneling process. 
In both parameter regions, the Landau-Zener adiabatic 
criteria in the double well can be easily satisfied by re- 
quiring ^A = 2ei J 2 along the horizontal loop in Fig. lc. 
Along the vertical loop, we choose ^ = 62A 2 to avoid 
the tunneling between two wells with different chemical 
potentials. Here e\ =0.1 <C 1, = 0.01 <C 1 are the 
adiabatic parameters. The energy gap is also numerically 
calculated and plotted in Fig. 2 for general parameters. 
We see the gap is always larger than \[2J for different 
A and U. When the lattice parameters are modulated 
along the elliptical loop in Fig. lc, there is no overall 
adiabaticity among double wells because atoms can now 
diffuse to a long distance if a long time period is used 
for the parameter modulation. Since the change of the 
parameters is slow, the cycle frequency of the periodic 
lattice modulation is much smaller than the energy gap. 

In the case of one atom per double well, the Hamilto- 




FIG. 5. (Color online) Plot of the diffusion of atoms after one 
cycle of the parameter modulation. A is fixed at 0. 5 varies 
from J to — J and then back to J. Solid line: initial density 
distribution. Dashed line: final density distribution after one 
cycle, (a) U = 5 J. (b) U = 20J. 



nian in each double well can be written as 



A 

-J 



-J 
-A 



(4) 



on the basis {|10),|01)}. The eigenenergies are E± = 
±a/A 2 + J 2 and the minimum energy gap is 2 J, therefore 
the above rate of the change of the lattice parameters still 
keeps the adiabaticity in each double well. 



III. QUANTUM TRANSPORT IN DOUBLE 
WELL OPTICAL LATTICES 

We numerically integrate the time-dependent 
Schrodinger equation with the Hamiltonian (pQ) us- 
ing the TEBD algorithm and calculate the density 
distribution of atoms at different lattice sites. In Fig. 3, 
we plot the number of atoms at lattice sites 15 and 16 
after one cycle of the parameter modulation described 
in Fig. lc. The lattice system contains total N = 20 
sites with 10 atoms uniformly distributed between sites 
5 and 14. Clearly the number of atoms transported to 
the neighboring double well (sites 15 and 16) depends 
on the interaction and the parameter modulation loop. 
In the weak interaction region (U < A), two atoms are 
transferred to sites (15, 16), which is the same as that for 
non- interacting bosons. In the case of non-interacting 
bosons, atoms follow the changes of the potential minima 
in the lattice system and it is easy to see that atoms 
would like to move one double well distance (two lattice 
sites) after one cycle of the parameter modulation. As 
the interaction increases, there is a sharp transition for 
the square loop, where the number of the transported 
atoms at sites 15 and 16 becomes one (Fig. 3a). In this 
region, the atoms behave the same as non-interacting 
fermions, and two atoms cannot occupy the same lattice 
site. The same physical picture applies to more cycles of 
the parameter modulation. In Fig. 4, we plot the density 
distribution of atoms after 10 cvcles of the parameter 
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FIG. 6. (Color online) Plot of the center of mass motion of the 
atoms in one cycle of the potential modulation. Initially six 
atoms are uniformly distributed in three double wells. Solid 
line: U = 5J, square loop. Dashed line: U = 5J, elliptical 
loop. Dotted line: U = 20J, square loop. Dashed dotted line: 
U = 20J, elliptical loop. 
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FIG. 7. (Color online) (a) Plot of the center of mass motion of 
the atoms in one cycle of the potential modulation. Initially 
six atoms are uniformly distributed in six double wells. Solid 
line: U = 5J, square loop. Dashed line: U = 5J, elliptical 
loop. Dotted line: U = 20J, square loop. Dashed dotted line: 
U — 20J, elliptical loop. 



modulation. For a small interaction strength, the atoms 
are transported to the target region that are 20 lattice 
sites away from the original region (Fig. 4a). For a large 
interaction strength, atoms move to both left and right 
directions (Fig. 4b), and the center of mass of the atoms 
does not change (i.e., no mass transport). 

When the parameter modulation is along the elliptical 
loop, there is a diffusion of the atoms along the lattices 
because the locally confined initial state is not the ground 
state of the Hamiltonian. In Fig. 3b, we plot of the num- 
ber of atoms at lattice sites 15 and 16 after one cycle of 
the parameter modulation along the elliptic loop. The 
main feature of Fig. 3a is still kept. In the weakly inter- 
acting region, the number of transported atoms is two, 
but the sharp transition is smoothed out and the num- 
ber of atoms at the large U is no longer one because of 
the diffusion of the atoms through the lattice. In Fig. 5, 
we plot the density distribution of atoms after one cy- 
cle of the parameter modulation with a fixed A = and 
varying 5 between 1 and —1. The diffusion of the atom 
density depends strongly on the rate of the change of 
the parameters. To suppress such uncontrolled diffusion, 
henceforth we mainly consider the square loop. 

The atom transport in the optical lattice can also be 
described using the center of the mass (COM) motion 
of the atoms. In Fig. 6, we plot the COM motion of 
atoms in the double well lattice after one cycle of the 
parameter modulation. We see in the weakly interacting 
region, the COM is shifted by two lattice site for both 
elliptical and square loops. In the strongly interacting 
region, the COM does not change for the square loop, 
indicating no mass transport, but varies for the elliptical 
loop, showing the asymmetry of the diffusion process in 
the parameter modulation process. 

In the above discussion, we consider that there are two 
atoms per double well (i.e., integer filling) in the initial 
state. Another interesting region is the half filling, that 
is, one atom per double well. In Fig. 7, we plot the 
COM motion of atoms after one cycle of the parame- 
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FIG. 8. (Color online) The transport of atoms in the presence 
of a harmonic trapping potential V% = Vo (i — x ) 2 , where x Q 
locates at the center of the lattice. 

ter modulation at the half filling. In both weakly and 
strongly interacting regions, the shift of the COM is two 
lattice sites, indicating that the non-interacting bosons 
and fermions behaves similarly at the half filling although 
their transport properties are completely different at the 
integer filling. The elliptical and square loops also yield 
similar results. 

In addition to the double well optical lattice potential, 
the atoms may also experience an overall harmonic trap- 
ping potential. In Fig. 8, we plot the transport of atoms 
in the presence of an harmonic trap. Clearly the har- 
monic trapping potential does not affect the transport 
even far away from the trap center. More interestingly, 
the initial right-moving COM motion turns to the left- 
moving motion after it reaches the maximum position 
that is determined by the harmonic trapping frequency. 
In a realistic experiment, this phenomenon corresponds 
to a dipole oscillation of atoms generated by the periodic 
modulation of the lattice parameters. The turnaround 
can be understood from the fact that when the COM 
reaches certain x, the chemical potential difference 2A 
cannot overcome the potential difference between two 
wells induced by the harmonic trap. Therefore atoms 
do not move at that cycle. Further modulation of the 
lattice parameters then provides a driving of atoms to 
the left-moving direction. 

Finally, we discuss how to observe the quantized atom 
transport in the double well optical lattices. In princi- 
ple, the transport dynamics can be observed using the 
single atom detection technology demonstrated recently 
|43l-l46|. Here we consider the signature of the quantum 
transport in more conventional experimental techniques: 
the momentum distribution in the time-of- flight image. 
The momentum distribution of the atoms in the optical 
lattices can be calculated as 

n k , t = 53<V(*)|oto J #(t))e i *< i -^. (5) 
The numerical results are shown in Fig. [9] for U = 5 J 



6 



(a) and U = 20 J (b). In the weak- interacting case, 
the atoms in one double well are perfectly transported 
to another double well and the final density distribution 
just have a global shift of the position, there we expect 
n>k,t — n>k,t+T after one period, as shown in Fig. [9^. In 
the strong-interacting case, the cold atoms move along 
two opposite directions, therefore rik,t 7^ n>k,t+T and the 
coherent peaks in the momentum distribution disappear. 
We find the momentum distribution is generally inde- 
pendent of the loops used in TEBD calculations. In ex- 
periments, the momentum distribution can be measured 
directly from the time-of- flight image, thus provides a di- 
rect experimental signature for the quantum transport. 



IV. A SEMICLASSICAL THEORY 



teracting bosonic atoms are equivalent to non-interacting 
fermionic atoms, therefore their transport properties are 
the same as the non-interacting electrons. For simplicity, 
we consider a periodic system to avoid the diffusion. 

With a periodic boundary condition, we can transform 
the Hamiltonian (pQ) without the interaction (U = 0) to 
the quasimomentum space 



H = r (q) ■ a 



(6) 



transformation 



d2j-l 

T,q eiq ' jaa qh Where 

and 



by using the Fourier 

^ e iq-(2j-l)a/2 aq ^ ^ 

T(q) = (-Jcos^,(5sin^, A), and the spin up 
and down correspond to the left and right sites of the 
double well. The Hamiltonian (|6j) has two energy bands 
a = zb with the dispersion 



Physically, the numerical results obtained in the above 
section at the half filling may be understood through the 
change of the Wannier function of bosons in a periodic 
lattice due to the modulation of the lattice parameters. 
The center of the Wannier function follows the poten- 
tial minimum of the double well, and can move from one 
double well to the neighboring one in one cycle. In this 
section, we present a semiclassical theory for the quan- 
tum transport of the non- interacting bosons in the double 
well lattice to understand the numerical results observed 
in the above section. On the other hand, the strongly in- 
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FIG. 9. Plot of the momentum distribution of atoms in the 
time of flight image, (a) U = 5J and (b) U — 20J. 
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The velocity of an atom in the bands satisfies the semi- 
classical equation of motion [47j 



Xa (q,t) = ~ Mat 



dq 



L qti 



(8) 



where fl^ t = -2lm(^^\^§fj is the a-th band Berry 
curvature in the momentum and time spaces, 



V2e a (e a - A) 



(Jcosf +i<5sinf ) 



(9) 



is the eigenwavefunction. The total particle transport 
along a close loop in the parameter space can be written 

as 



c T 



j (t) dt 



-E 

a=± 



7r I a 



7v /a 



dt I f(e a - n) x a (q, t) dq, 
(10) 

where j (t) is the atom number current, / (e a — fi) 
is the Fermi-Dirac distribution / (e a — fx) = 
1/ [exp ((e a — ji) /ksT) + 1] for fermionic atoms 
and the Bose-Einstein distribution f (s a ) — 
1/ [ex.-p((e a — /i) /ksT) — 1] for bosonic atoms. For 
the non-interacting bosonic atoms, all atoms occupy 
the lowest energy state and the chemical potential \i 
loses its meaning and total atom transport at the zero 
temperature becomes 



c T 



^ j> dtX- (<2min,£) , 



(11) 



where only the lowest energy band has a contribution, 
<2min is the quasimomentum at the energy minimum, r] 
is the number of atoms in each double well. The energy 
minimum of the lowest band E- locates at q m { n = ex- 
cept for the parameter S = ±1, where the band is flat and 
does not depend on q. Note that 5 = ±1 corresponds to 
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isolated double wells in the lattice (i.e., no inter- well tun- 
neling), therefore the initial state of the system can be 
chosen with q = 0. 

In the case of non-interacting bosonic atoms with an 
initial q = 0, the first term in x_, ^-\ q= o = 0, and the 
second term 



ot 4 dt ( J2 -f A 2 - AVJ 2 + A 2 ) ' V ; 

Therefore the total atom transport along the square loop 
in Fig. lc is 

ct = — — j> dtVt^ — r] (13) 

when A ^> J. When there are rj atoms in each double 
well, there are n atoms pumped out the system, agreeing 
with the numerical results presented in Figs. 6 and 7. 

In the case of strongly interacting bosons, which is 
equivalent to non- interacting fermions, the atoms grad- 
ually occupy different q states as the number of filling 
increases. In the case of the half filling, the lowest energy 
band S- is fully occupied and the chemical potential /i 
lies at the gap between two bands. The total transported 
atoms can be shown to be ct = 1 using Eq. ([TO]) [5j, 
which agrees with the numerical results obtained in Fig 
7. In the case of the integer filling (two atoms per double 
well), both bands are fully occupied and ct — because 
the Berry curvatures = —Q~ t for two bands and their 
contributions to ct cancel with each other. This result 
agrees with the numerical results presented in Fig. 6. 

Finally, we comment on the validity of using the pe- 
riodic boundary condition in the semiclassical theory for 
the explanation of the numerical results observed in Sec. 
III. Eqs. ([6ll7([TT]) are introduced by assuming a periodic 
boundary condition for a single atom and the fact that 
the atom occupies the Bloch eigenstate. In general, the 
periodic system is not equivalent to our system where 
the atoms are initially localized. However, because the 
initial state we choose is the product of the state in each 
isolated double well (no inter- well tunneling), and the 



atoms at different double wells do not affect each other 
(no interaction), the periodic system can be taken as a 
multiply copies of our local system. Therefore the ef- 
fect of the quantum transport and the essential physics 
should be the same for the periodic system and our local 
system, which justify the agreement between the semi- 
classical theory and the numerical results. The difference 
between these two systems is that, in a periodic system, 
the change of the local atom density after a period cannot 
be observed because it is periodic. Instead, what can be 
observed is the current, whose integration in one period 
gives the number of transported atoms. While in our lo- 
cal system, the transport of atoms is directly reflected 
on the density variation because the neighboring lattice 
sites are not occupied initially. 

V. CONCLUSION 

In conclusion, in this paper we study the quantum 
transport of ultra-cold bosonic atoms in double well op- 
tical lattices where the lattice parameters can be peri- 
odically modulated. The transport of atoms depends 
strongly on the atom filling and the interaction. In the 
strongly interacting region, the bosonic atoms behave 
similarly as fermions with quantized transport at the 
half filling and no transport at the integer filling. In 
the weak interacting region, the quantized transport are 
robust and the number of the pumped atoms is propor- 
tional to the atom filling per double well. Signature of the 
quantum transport of atoms in the momentum distribu- 
tion is obtained. In addition to the numerical simulation 
of the transport dynamics of atoms in the double well 
optical lattice using the TEBD algorithm, we develop a 
semiclassical model to explain the numerical results in 
the non-interacting and strongly interacting limits. Our 
scheme for the quantized atom transport does not in- 
volve accurate control of the lattice parameters and their 
time dependence, thus may provide a robust way for dis- 
tributing atoms in optical lattices, which are critically 
important for quantum computation in optical lattices 
as well as the atomtronics applications. 
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